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Abstract. Differences in cloud droplet effective radius and cloud droplet number concentration (CDNC) estimates inferred 
from the Aqua MODIS Collections 5.1 and 6 cloud products (MYD06) are examined for warm clouds over global oceans for 
the year 2008. Individual pixel level retrievals for both collections are aggregated to 
1° x 1° and compared globally and regionally for the three main spectral channel pairs used for MODIS cloud optical 
property retrievals. Comparisons between both collections are performed for cases in which all three effective radii retrievals 
are classified by the MODIS Cloud Product as valid. The contribution to the observed differences of several key MYD06 
Collection 6 algorithm updates are also explored, with a focus on changes to the surface reflectance model, assumed solar 
irradiance, above cloud emission, cloud top pressure, and pixel registration. Global results show a neutral to positive (> 50 
em”) change for C6-derived CDNC relative to C5.1 for the 1.6 um and 2.1 pm channel retrievals, corresponding to a neutral 
to -2 um difference in droplet effective radius. For 3.7 um retrievals, CDNC results show a negative change in the tropics, 
with differences transitioning toward positive values with increasing latitude spanning -25 to +50 cm” related to a +2.5 to -1 
um transition in effective radius. Cloud optical thickness differences were small relative to effective radius, and found to not 
significantly impact CDNC estimates. Regionally, the magnitude and behavior of the annual CDNC cycle are compared for 
each effective radius retrieval. Results from this study indicate significant intercollection differences in aggregated values of 
effective radius due to changes to the pre-computed retrieval lookup tables for ocean scenes, changes to retrieved cloud top 
pressure, solar irradiance, or above cloud thermal emission, depending upon spectral channel. The observed differences 
between collections may have implications for existing MODIS derived climatologies and validation studies of effective 


radius and CDNC. 


1 Introduction 

Marine stratocumuli play a large role in the modulation of the Earth’s radiative balance and hydrological cycle, 
owing to their persistence and large areal extent. They are maintained by moisture flux from the sea surface, longwave 
cooling, and dry air entrainment at the cloud top. Confined to the planetary boundary layer, most often in subtropical 


subsidence areas, they are susceptible to the influences of anthropogenic aerosols. Aerosol perturbation manifests changes 
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though several mechanisms. Through the first indirect aerosol effect, anthropogenic influences potentially alter the degree 
of reflection of shortwave radiation by distributing a given quantity of cloud water over a larger number of droplets 
compared to an unperturbed cloud by supplying additional cloud condensation nuclei (CCN) (Twomey, 1974). This increase 
in cloud albedo provides a local cooling effect as the cloud top temperature is near the sea surface temperature, producing a 
negligible longwave radiative effect. With a decrease in droplet effective radius from the addition of CCN, autoconversion 
efficiency can decrease reducing precipitation efficiency (Albrecht, 1989) and potentially increasing the overall lifetime of 
clouds, augmenting the first indirect effect (Lohmann and Feichter, 2005). This naturally leads to a need for estimates of 
cloud droplet number concentration (CDNC) from observation to better understand the natural background and 
anthropogenic contributions to cloud droplet number concentration. 

Satellites offer an opportunity to address this need through their global observational area spanning decades. The 
National Aeronautics and Space Administration’s (NASA) Moderate Resolution Imaging Spectroradiometer (MODIS) is an 
instrument suited for this task. Flown aboard two Earth-observing satellites (Terra and Aqua), MODIS provides 
observations in the appropriate spectral bands to infer cloud optical thickness and droplet effective radius through the bi- 
spectral method of Nakajima and King (1990). From the retrieved cloud optical thickness and droplet effective radius, cloud 
droplet number concentration can be estimated (Brenguier et al., 2000). NASA provides retrievals of the cloud optical 
properties through the MODIS Cloud Product (MOD06 and MYD06 for Terra and Aqua, respectively); note the MOD06 
and MYD06 products are produced by the same algorithm (Platnick et al., 2015; Platnick et al., 2017). As passive imager 
cloud remote sensing science evolves, the MODIS Cloud Product undergoes periodic updates to the cloud screening and 
retrieval algorithms, which are collectively implemented and reprocessed as “collections.” Collection 6 is the most recent 
release of the Cloud Product and includes several changes that propagate through to estimates of cloud microphysical 
properties relative to its predecessor, Collection 5.1 (Platnick et al., 2015). Since there is an existing body of research based 
upon Collection 5.1 of effective radius (e.g. Zhang and Platnick, 2011; Painemal and Zuidema, 2011) and CDNC (e.g. Zeng 
et al., 2014; Ahmad et al., 2013) there is a need to understand how these retrievals differ between collections. In this study, 
we investigate and highlight the differences in retrieved cloud microphysical properties manifested both in effective radius 
and CDNC between both collections using one year of Aqua MODIS (MYD06) observations of warm clouds over global 


oceans in an attempt to understand the differences in CDNC from the respective products. 


2 Data and Methods 

The MODIS Collection 5.1 (C5.1) Cloud Product provides estimates of cloud optical thickness (z) and droplet 
effective radius (r,) in addition to cloud screening, quality control flags, and ancillary data relevant to the computation of 
cloud droplet number concentration. For most parameters derived from visible and near/shortwave infrared channels, the 
spatial resolution is | x 1 km at nadir. Some parameters, such as those derived from thermal infrared channels, are provided 


ata 5x5 kmresolution. The C5.1 Cloud Product applies a Clear Sky Restoral (CSR) algorithm that identifies pixels that are 
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expected to be only partially cloudy, and excludes them from the cloud optical properties retrieval, essentially limiting C5.1 
data to cloudy scenes only. In addition, the C5.1 r. retrievals using the 1.6 and 3.7 um channels are given as a difference 
relative to the 2.1 um 7, value, and are therefore dependent upon the success of the 2.1 um retrieval. 

The Collection 6 (C6) Cloud Product offers improvements in retrievals of several relevant parameters for CDNC 
estimation. Changes to the optical and microphysical retrievals for C6 include improvements to the forward radiative 
transfer model used to create the pre-calculated cloud retrieval look-up tables, re-registration of the visible/near IR focal 
planes for MODIS Aqua, and explicit reporting of all three spectral channel pair effective radius retrievals (Platnick et al., 
2015); cloud top property retrieval changes for C6 include additional cloud top temperature and pressure retrievals at 1 x | 
km spatial resolution and improvements to the characterization of low-level clouds over water surfaces, among others 
(Platnick et al., 2015). With respect to the assessment of pixel quality, the CSR algorithm is still applied and its results are 
reported in the cloud product, but unlike C5.1, cloud optical property retrievals are attempted on pixels identified as partly 
cloudy and, if successful, are written to the cloud product separate from the heritage cloudy retrievals (Platnick et al., 2015). 
Note that pixels identified by the CSR algorithm as being not cloudy (e.g., false cloudy sun glint or thick aerosols such as 
smoke or dust) remain excluded from the cloud optical property retrievals. One year of MODIS Aqua C5.1 and C6 data, 


corresponding to calendar year 2008, is used in this investigation. 


2.1 Calculating CDNC 

Relating cloud droplet number concentration to retrievals of optical thickness and droplet effective radius requires 
an estimate of cloud liquid water content. To attain CDNC, which is not a parameter in the MODIS cloud product, the 
vertical structure is assumed to follow the so-called adiabatically-stratified cloud model (Brenguier et al., 2000). This model 
accounts for a linear increase in cloud liquid water content as a saturated air parcel rises through a cloud’s vertical extent 
while undergoing droplet growth as the parcel cools moist-adiabatically. As marine boundary layer clouds are typically 
shallow, the temperature is nearly constant, so it follows that the amount of condensate formed through adiabatic ascent can 
be assumed to be constant. The following relationships, adapted from Bennartz (2007), are presented as a summary of how 
CDNC relates to the adiabatic cloud model. The adiabatic liquid water content w4p at any point 4 above the cloud base can 


be expressed as: 


W,(h,T)=c¢,(T )h (1) 


where c,, is the mass of condensate formed per unit volume for each meter of ascent and T is the cloud top temperature. The 
Cy parameter is derived from the Clausius-Clapeyron relationship, and is primarily a function of temperature and, to a lesser 


degree, pressure. For this study, an assumed cloud-top pressure of 850 hPa is used, which corresponds to a geopotential 
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height of approximately 1500 m. For marine stratocumulus regions such as the Southeast Pacific, this height is near the 
mean top of the boundary layer (von Engeln and Teixeira, 2013) around Aqua’s 13:30 local equatorial crossing time (Ho et 
al., 2015). 

Relating the liquid water content profile to the cloud optical properties, Brenguier et al. (2000) express the adiabatic liquid 


water path as: 


(2) 


where the effective radius is the value at cloud top. Accounting for cloud-top entrainment of dry air, observed and modeled 
liquid water paths are typically near 80% of the purely adiabatic value (Duynkerke et al., 2004; Pawlowska and Brenguier, 
2003). The c,, parameter in this study is therefore scaled to approximately 80% of the purely adiabatic value, in order to 
avoid overestimation of the condensate observed by MODIS. 


With the cloud optical parameters, Bennartz (2007) expresses CDNC as: 


T° -s2| 3 i ess - ne 
N= —-[2W0] [20] =| C, (3) 


where k encapsulates the skewness and dispersion of the cloud droplet size distribution, Q is the scattering efficiency, and p; 


3 
r 
is the density of liquid water. Analytically, k = = , where 7, is the volume mean droplet radius. For marine boundary 


a 


layer clouds k=0.8 is assumed to be a representative value for this study (Brenguier et al., 2011; Martin et al., 1994), 
although Painemal and Zuidema (2011) suggest that A=0.88 may be more representative for droplets near the cloud top, 
resulting in a narrower cloud droplet size distribution. As the size parameters for cloud droplets approach the geometric 


optics limit, the asymptotic value of 2 is assumed for Q. 


3 Results 
3.1 Global-Scale Common Pixel Comparison 

Common pixel scenes are those in which all three effective radius retrievals are valid for both collections with 
otherwise consistent selection criteria. A common pixel selection allows for an objective comparison of effective radius and 
CDNC estimates between both collections, free of the influence of pixel population differences due to, for example, 


reclassification of cloud phase due to changes to the cloud thermodynamic phase algorithm (Marchant et al., 2016), or where 
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one or more effective radius retrievals have failed between collections (Cho et al., 2015). The global distribution of the total 
MODIS Aqua common pixel count for calendar year 2008 is shown in Fig. 1. 

Fig. 2 shows mean annual C6 values and C6-C5.1 differences in MYD06 cloud droplet effective radius for common 
pixel scenes, aggregated globally to a 1° x 1° equal angle grid for calendar year 2008, from the 1.6, 2.1 and 3.7 um channels. 
For r.;.6, C6 retrievals are smaller than C5.1 over most of the global oceans (differences generally less than | zm), with the 
largest differences (greater than 2 um) occurring near the Baja California peninsula. In equatorial regions of the Indian and 
West Pacific Basins, as well as over sea ice regions in the northern high latitudes above 60° N, however, 7.1.5 differences are 
largely neutral to slightly positive. C6 retrievals of 7.2; are predominantly smaller than C5.1 (differences less than | um), 
with the greatest disagreement between collections occurring in the tropics. For 7.37, there is a latitudinal gradient in the 
retrieval differences, with positive intercollection differences in the tropics (C6 > C5.1 by more than 2 ym) transitioning to 
slightly negative values (C6 < CS5.1 by less than 2 um) at higher latitudes. 

The relationship between the intercollection shifts in r, do not necessarily imply that the cloud liquid water path 
profile is steepened in C6 by the decrease in r.;6and 7-2, and the increase of 7.37 at low latitudes. Contributing effective 
radii in this study are not required to conform to an adiabatically stratified profile and are likely subject to retrieval artifacts, 
complicating the determination of the liquid water content profile. This is especially the case at lower latitudes where 
convective clouds are more dominant and inhomogeneous relative to stratiform clouds. 

Patterns similar to those of effective radius are mirrored in the respective CDNC C6 average and intercollection 
difference plots of Figure 3. As shown in Fig. 3d-e, for the N;.5 and N2, estimates, the differences are generally neutral to 
positive, i.e., C6 > C5.1, with values over +50 cm™ occurring along coastal boundaries in the subtropical subsidence zones. 
Fig. 3f shows N37 from C6 is smaller than C5.1 in the tropics and generally larger at latitudes higher that 45°. Note that of 
the three spectral channels used for 7, retrievals, the sensitivity of the 3.7m channel is weighted closest to cloud top 
(Platnick, 2000), and r.37 is therefore considered to be the most appropriate effective radius for estimating CDNC based 
upon the assumption of the adiabatic model that the effective radius is the cloud top value. 

Cloud optical thickness, which is another necessary retrieved parameter for estimation of CDNC, is neutral to 
slightly higher in aggregated retrievals for Collection 6 (not shown). For regions in which stratocumuli are often present, 
generally in the subtropical subsidence zones, C6 cloud optical thickness is within 0.25 of C5.1. This translates to a relative 
overestimation by C6 of roughly 2%. Considering the relatively weak sensitivity of CDNC to changes in cloud optical 


thickness relative to droplet effective radius, this results in an increase of C6 CDNC by less than 1%. 


3.2 Algorithmic and Calibration Factors in Intercollection Differences 
Due to the strong sensitivity to effective radius (N «r;°” ) in Equation 3, larger effective radii translate into lower 


relative CDNC (and vice versa), thus the r, differences shown in Fig. 2 largely explain the CDNC differences observed in 
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Fig. 3. As previously noted, the use of common pixels in the analysis in Section 3.1 removes the influence of pixel 
population differences between the two collections. The large effective radius differences between C5.1 and C6 observed in 
Fig. 2 are therefore the result of either changes to the MYD06 optical property retrieval algorithm itself or changes to 
upstream or ancillary products such as the Level 1b calibrated radiometric data or the MYD06 cloud top property retrieval 
algorithm. 

For C6, changes to the liquid phase cloud optical property retrieval algorithm were, for the case of 7.7.6 and r..21, 
limited primarily to the pre-computed retrieval look-up tables (LUTs). The change having the most significant impact on the 
re retrievals is the use of new LUTs for retrievals over ocean, which incorporate an ocean surface bidirectional reflectance 
model directly into the LUT that uses the Cox-Munk wind speed and direction-dependent wave-slope distribution (Cox and 
Munk, 1954b, a). Previously in C5.1 the same LUT was used over land and ocean scenes, with the ocean surface assumed to 
be a Lambertian reflector with 5% albedo (Platnick et al., 2017). Changes to input datasets, such as Aqua Level 1b re- 
registration of the visible/near-IR focal plane and updates to the cloud top property retrievals, are also found to significantly 
impact the r, retrievals. Figs. 5, 6, and 7 show a single granule case study illustrating the impacts of these changes. An Aqua 
MODIS granule of open and closed cell marine stratocumuli acquired off the coast of Baja California on 9 June 2014 (2155 
UTC) was selected. The true color image, 1km CTP, 7, and 7,2; are shown in Fig. 5a-d. A region of sunglint is visible 
slightly west of nadir below 15°N. Also, a glory is present east of nadir above 23°N, which is evident in the difference plots 
in Fig 6. The C5.1 to C6 re.7.6, 22.1, and re.3,7changes, defined here as the difference r,.,(C6) — 7.,(C5.1), are shown in Fig. 
6a-c. Note that purple colors arise in regions where adjacent pixels have positive (red) and negative (blue) differences in 
effective radii. Subsequent rows in Fig. 6 show the individual impacts of the C6 changes listed above. 

The impacts of the new LUT for ocean scenes on r..;.6, %¢2.;, and 7.3.7 are shown in Fig. 6d-f. This change primarily 
impacts r, retrievals for optically thinner clouds (z roughly < 2 to 3) and those over sun glint, lowering r, by up to 1 micron. 
The strength of this effect is dependent on the orthogonality of the solution space, and thus the impacts are largest for re.1.5 
and smallest for r.37. These r, decreases due to LUT changes can be offset or enhanced, depending on the location of the 
pixel relative to cloud illumination or shadowing, by the effects of the Aqua Level 1b re-registration of the visible/near-IR 
focal plane, shown in Fig. 6g-i that are again strongest on r,;6 due to enhanced LUT non-orthogonality (7.3.7 is largely 
unaffected). The largest contribution to the 7... ; changes observed in Fig. 2 result from the updates to the cloud top pressure 
(CTP) retrievals, which are discussed in detail in Baum et al. (2012). The impacts on each channel are shown in Fig. 6j-1. For 
the liquid phase clouds in this granule, C6 CTP is largely higher than in C5.1 with increases of 100 hPa or more (not shown). 
Higher CTP, i.e., lower cloud top height, results in larger above-cloud atmospheric gaseous absorption corrections, brighter 
top of cloud reflectance, and thus smaller r, from all three spectral channels, by up to 1 micron or more. The global yearly 
mean CTP change for 2008 is shown in Fig. 4a, again illustrating that C6 CTP is largely higher than in C5.1, particularly 


over the marine stratocumulus regions where CTP have been increased to better place the height of the capping inversions 
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(Baum et al., 2012). In contrast, in regions where the liquid cloud occurrence is low (e.g., the tropics; see Fig. 1) C6 CTP 
slightly decreases. 

While the above C6 changes largely explain the intercollection differences for 7.5 and 7.2, algorithm updates for 
the 7,37 retrieval were more extensive, involving fixes to known shortcomings in the C5.1 algorithm, and their net 
contribution to the dissimilar difference patterns shown in Fig. 2 are unclear. Retrieved 7.37 decreases resulting from the 
CTP changes are enhanced by a correction to the assumed band-averaged solar irradiance (Fo) at 3.7um. For C6, Fo is 
smaller than for C5.1 and earlier versions (10.9 Wm-2um-1 for C6 versus 11.7 Wm-2um-1 for C5.1), resulting in a decrease 
of r.3.7 by 1 micron or more, as shown in Fig. 7a. However, above-cloud atmospheric emission, previously ignored in C5.1, 
is now accounted for in C6 and yields sizably larger r.37, on the order of 1 micron or more, as shown in Fig. 7b. The 
cumulative effects of the C6 changes discussed here appear to explain the observed granule-level intercollection 7,37 
differences (Fig. 6c), which for the case study granule are strongly positive (C6 larger) or negative (C6 smaller) and are 
seemingly dependent on view angle. Nevertheless, it is difficult to determine how the granule-level 7,37 differences translate 
to the global aggregated differences in Fig. 2. In addition, given the extensive C6 changes and the known shortcomings of 
the C5.1 7,37 retrievals it is of limited benefit to further assess their intercollection differences. 

A final consideration regarding the intercollection CDNC differences shown in Fig. 3 is the fact that CDNC 
estimates from the retrieved cloud optical properties also depend on the retrieved cloud-top temperature (CTT) through the 
Cy parameter. As discussed above and in Section 2, a number of algorithm improvements were implemented in the C6 
MYD06 cloud-top property product, including cloud-top retrievals at 1 km resolution in addition to the 5km cloud-top 
retrievals that have heritage to C5.1 and earlier collections (Baum et al., 2012). For the C6 CDNC estimates in this study, 
the 1 km CTT is used. To assess the impact of differences in CTT on CDNC, the CTT values from C5.1 at 5 km and C6 at 1 
km were aggregated to 1° x 1° for the common pixel population, shown in Fig. 4b. The differences in CTT reveal a 
relatively lower temperature (3K) from C6 in subsidence regions at low latitudes, neutral to slightly positive differences at 
higher latitudes, and largely positive differences in the intertropical convergence zone (ITCZ). Although the CTP increases 
in subsidence regions between collections, CTP placement is independent of CTT changes between collections, which 
creates seemingly contradictory situation when viewed together. To determine how the disagreement in CTT translates to 
CDNC apart from the cloud-top retrieval impacts to r,, the c,, was also aggregated to 1° x 1°. The differences in c,,5 and 
Cys translate to a 0.5% relative decrease in globally averaged N for C6, provided all other parameters are held constant. 
The largest observed difference is in the subtropical subsidence regions, where the decrease is near 4%. Therefore, c,, has 
only a marginal impact on N for any of the three absorption channels, a result that provides further evidence in favor of 


attributing intercollection N differences largely to changes in the effective radius retrievals. 
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3.3 Regional Common Pixel Comparison 

Four maritime regions were selected for additional examination primarily for their proximity to sources of 
anthropogenic aerosols and generally high number of liquid cloud observations within their domains: Southern Africa 
(SAF), South America (SAM), North America (NAM), and Southeast Asia (SEA). The boundaries of these regions are 
shown as white outlines in Fig. 1. 

The annual CDNC cycles from C5.1 and C6 are given in Fig. 8a for the Southern Africa region. Similar to the 
global results, the 1.6 and 2.1 um CDNC from C6 are systematically higher (+10 em>), while the C6 3.7 um values vary 
from being lower (< 10 cm™) to near agreement with C5.1. During the Southern African dry biomass burning season 
between June and October (Roberts et al., 2009), the 1.6 um C6 CDNC presents the greatest differences between collections, 
exceeding 30 em” with a maximum CDNC of over 180 cm”. Interestingly, there is no discernable difference at 3.7 um for 
either collection over the biomass burning season, however it is quite muted relative to N;.5; peak values during the season 
are roughly in agreement with Austral summer values between 90-100 cm”. For the region, absorbing aerosols from 
biomass burning often overlie the cloud deck, decreasing the reflected radiation observed by MODIS. With increasing 
aerosol optical depth, this can result in decreased retrieved cloud optical thickness and a general decrease in retrieved 
effective radius (Haywood et al., 2004). The sensitivity to overlying aerosols on the retrieved quantities depends on the 
degree of orthogonality of the retrieval solution space. Since CDNC estimates are significantly more sensitive to r, changes 
than optical thickness, the decrease in r, may result in an overestimation of CDNC if overlying absorbing aerosols are 
present. Results from the SAM region in Fig. 8b show a stratification and general response between collections, which is 
similar to SAF. For this region, a generally consistent aerosol emission rate gives rise to little variability (Huneeus et al., 
2006) in CDNC, with annual cycle amplitudes limited to a range of approximately 20 cm” across all channels for both 
collections. For NAM, the common pixel comparison also yields similar results as the previous regions as shown in Fig. 8c. 
For this region, both the 1.6 and 2.1 um channels also exhibit a 20-30 cm® difference for Collection 6 retrievals while 3.7 
um is in near agreement with Collection 5.1. All three channels have a uniform response to increases in CDNC coincident 
with a relatively early wildfire season that occurred during the spring and summer of 2008 in California (Brioude et al., 
2009). 

In contrast to the other regions in the study, the SEA domain yields a different stratification of CDNC shown in Fig. 
8d, where N37 > N).6 > N2, for both collections with general agreement between N37 retrievals. Disagreements between 
collections are less than 15 cm® for each channel. The different stratification may be a result of differing cloud dynamics in 
this region due to a weaker inversion strength, resulting in greater entrainment compared to the other regions which lie in 
subsidence zones. There is a strong seasonal response to aerosols visible in each channel coincident with the advection of 


aerosols by the Asian winter monsoon (Bennartz et al., 2011). 
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The adiabatic model assumes the retrieved effective radius is the near cloud-top value. For an ideal adiabatically 
stratified cloud 7.37> re21> e216, implying that Nj.6 > N2; > N37, as each effective radius retrieval has a varying vertical 
penetration depth. None of the focus regions in this study exhibit this behavior. In addition to entrainment processes, 
retrieved droplet effective radius and therefore estimated CDNC can be influenced by drizzle processes (Suzuki et al., 2010) 
or 3-D radiative transfer effects and cloud field inhomogeneity at the sub-pixel level that violate the plane-parallel cloud 
assumption used in the optical properties retrieval (Zhang and Platnick, 2011). Since values in this study are aggregated to | 
x | degree and averaged monthly with no screening criteria with respect to vertical stratification, it is likely these processes 


are contributing to the non-adiabatic stratification observed. 


4. Conclusions 


It has previously been shown, under the condition of adiabatic clouds, that estimates of the CDNC of warm marine 
liquid-phase clouds can be derived from passive satellite remote sensing observations. To this end, observations from 
MODIS on the Terra and Aqua satellites has seen wide use. The operational MODIS Cloud Product (MYD06 for Aqua) 
provides the cloud property retrievals necessary for computing CDNC estimates, namely cloud-top temperature, cloud 
optical thickness, and effective droplet radius (7,). Collection 6, the most recent release of MYD06, includes numerous 
updates relative to its predecessor, Collection 5.1, that can propagate through to estimates of CDNC. Using one year (2008) 
of global MODIS Aqua observations, intercollection differences are investigated for CDNC derived from three independent 
spectral r, retrievals, namely from the 1.6, 2.1, and 3.7um channels. For the pixel population having successful r, retrievals 
from all three spectral channels in both collections, C6 rz; and re) retrievals are smaller than those of C5.1, with the 
greatest differences generally found for r.;. These intercollection r, differences result in relatively larger estimates of C6 
CDNC for both channels. Nevertheless, CDNC from these two spectral channels offer similar annual cycles for both CS.1 
and C6, regardless of the geographic region. N37 differences in stratocumulus regions are generally smaller and, unlike N;5 
and N>,, the signs of the differences are rarely uniform throughout the annual cycle. Moreover, the global intercollection 
differences in r.37 and N37 exhibit quite different behavior than those derived from the 1.6 and 2.1m channels. These 
differences are inherently attributable to the more wholesale C6 changes that addressed known shortcomings in the C5.1 r..3.7 
retrieval. But because the individual effects of the C6 changes can be quite large, are often of opposite sign, and have 
angular and potentially other unknown dependencies, it remains unclear how granule-level r..3.7 differences translate to the 
global aggregated differences shown here. Additional research is necessary in order to quantify the contribution of these 
dependencies to the observed intercollection differences of r..37 and N37. Furthermore, it is recommended that quantitative 


use of the C5.1 r,.3.7 retrieval be avoided given its known shortcomings. 
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Among the updates for C6 is the independent reporting of r, for each wavelength, rather than as differences with 
respect to 7,2, aS well as the inclusion of retrievals of the so-called partly cloudy pixels previously discarded in C5.1. These 
changes offer additional options in the analysis of CDNC from MODIS. Several permutations in data screening scenarios 
beyond those used in this study are now available which were not possible in C5.1, and will likely be explored in future 
studies. Caution will be warranted, when interpreting the results of these and other future studies of r. and CDNC inter- 
comparisons which may be subject to different screening criteria available in C6 which allows for the inclusion of more 
potentially problematic retrievals. The differences observed here are for overcast scenes only for pixels in which all three 
effective radius retrievals are valid for liquid phase clouds. Different screening criteria may exacerbate the effects of 


aggregated retrieval artifacts, which have not been considered in this study. 
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Figure 1: Annual total of Aqua MODIS liquid phase cloud observations where all three effective radii retrievals are 
valid per 1 x 1 degree grid box. Focus regions are overlain. 
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Figure 2: Mean annual Collection 6 cloud droplet effective radius for calendar year 2008 for the a) 1.6 pm, b) 2.1 pm, c) 3.7 pm 


chagnels and C6 — C5.1 effective radius differences for the d) 1.6 1m, e) 2.1 pm and f) 3.7 wm channels. 
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Figure 3: Figure 2: Mean annual Collection 6 cloud droplet number concentration for year 2008 for the a) 1.6 1m, b) 2.1 pm, c) 3.7 


pum channels and C6 — C5.1 CDNC differences for the d) 1.6 pm, e) 2.1 pm and f) 3.7 pm channels.. 
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Figure 4: Cloud top pressure (a) and temperature (b) differences between aggregated 1 km resolution C6 and 5 km resolution C5.1 
for the common pixel population. 


16 


Cloud Top Pressure (mb) 


200 400 600 800 1000 
. 2 hd tJ dl 4 ie . 


re 
an - - i” o- - - 


Ice Phase COT Liquid Phase COT Ice Phase CER (um) Liquid Phase CER (em) 
i RS EEE 


ot 1 10 100 G1 1 10 msEBESSHSEHEHE SES PSA S 


Figure 5. (a) True color image, and C6 (b) 1km CTP, (c) cloud optical thickness, and (d) 2.1m effective radius for a 


MODIS Aqua granule acquired off the coast of Baja California on 9 June 2014 (2155 UTC). 
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Figure 6. C6-C5.1 r, retrieval differences for the granule shown in Fig. 5 for r.7.6, re2.1, and r,,37 effective radii a-c) in 


total and contributions from d-f) lookup table differences, g-i) pixel reaggregation, and j-l) cloud top pressure. 
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Figure 7. C6 — C5.1 r,.3.7 retrieval differences due to (a) the change of the assumed solar irradiance (Fo) at 3.7um and 


(b) correctly accounting for above-cloud atmospheric emission previously ignored in C5.1. 
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Figure 8: Annual CDNC cycle plots for the a) Southern African, b) South American, c) North American and d) Southeast Asian 
regions. 
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